function p = parameters(~)

% elasticity of substitution between tasks, from Humlum (2021)
p.lambda = 0.5;

% parameters of workers' productivity
p.m = 1; % m>0 for wages increasing in skill
p.n = 0;
p.a = 2; % a>0 for monotone comparative advantage across workers
p.A = 1; % scale parameter

% parameters of capital productivity
p.b = -20; % b<0 for (quasi)concavity of psi_{k,x}/psi_{s,x} in x
p.ak = 8;

% labor supply
% choose "baseline" for initial equilibrium (normalized to 1)
% choose "new" to change labor supply to the function labornew.m
p.labor = "baseline";
%p.labor = "new";

% capital price decline from 1980 to 2016 (used for computing calibration
% targets with targetfun.m)
% choose which price index to use
p.qindex = "baseline"; % non-residential equipment and software, quality-adjusted (baseline)
%p.qindex = "eqqa" % non-residential equipment, quality-adjusted
%p.qindex = "eq" % non-residential equipment
%p.qindex = "sw" % software

% parameters for solution algorithms
p.contol = 10^-10; % constraint tolerance for output maximization
p.stol = 10^-10; % tolerance for stilde: if stilde is within p.stol of 0 (or 1), treat it is 0 (or 1)
p.qtol = 10^-5; % tolerance for q: if q is within p.qtol of qinfinity, assume unbounded output
p.optionsfzero = optimset('TolX',eps); % default: TolX=eps, use TolX=10^-7 for speed

end